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We investigate theoretically and experimentally a nondestructive interferometric measurement of 
the state population of an ensemble of laser cooled and trapped atoms. This study is a step towards 
generation of (pseudo-) spin squeezing of cold atoms targeted at the improvement of the Caesium 
clock performance beyond the limit set by the quantum projection noise of atoms. We calculate 
the phase shift and the quantum noise of a near resonant optical probe pulse propagating through 
a cloud of cold 133 Cs atoms. We analyze the figure of merit for a quantum non-demolition (QND) 
measurement of the collective pseudo-spin and show that it can be expressed simply as a product 
of the ensemble optical density and the pulse integrated rate of the spontaneous emission caused by 
the off-resonant probe light. Based on this, we propose a protocol for the sequence of operations 
required to generate and utilize spin squeezing for the improved atomic clock performance via a QND 
measurement on the probe light. In the experimental part we demonstrate that the interferometric 
measurement of the atomic population can reach the sensitivity of the order of s/N a t in a cloud 
of N at cold atoms, which is an important benchmark towards the experimental realisation of the 
theoretically analyzed protocol. 
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I. INTRODUCTION 

The quantum features of a collective atomic magneti- 
zation were experimentally addressed for the first time 
in the study of quantum noise of a collective spin done 
by Alexandrov and Zapasskii [| . The relevance of this 
work was later on accentuated by the first observation of 
a quantum noise limited magnetization, i.e. linear depen- 
dence of the atomic variance on the number of atoms, by 
S0rensen et. al. [2|. This quantum limit, called projec- 
tion noise limit, has been reached in state-of-the-art Cs 
atomic clocks [3j and, by now, is the limitation towards 
the improvement of the clock precision. However, the 
projection noise would not be the limiting factor for the 
clock precision if one was able to increase the number of 
atoms used in the clock operation, but this has not been 
possible due to large collisional shift in cold Cs samples. 

Nevertheless, it is possible to overcome the projection 
noise limit, as it has been demonstrated with the gen- 
eration of entangled and squeezed states of two ions Q , 
and the creation of spin squeezed states in a cloud of 
cold excited Cs atoms 5J. As a protocol for the gen- 
eration of spin squeezed states, the use of a quantum 
nondemolition (QND) measurement has been proposed 
and implemented in a vapor cell in [6j, and the same 
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kind of interaction was first proposed in |y, |]j as a mean 
for improving the clock performance. Another approach 
towards performing atomic spectroscopy below the stan- 
dard quantum limit using Bose-Einstein condensates has 
been suggested in 

In this paper we propose a sequence of spin rotations 
and QND measurements which should allow to overcome 
the projection noise limit for a Cs clock. The method 
involves generation of the coherent superposition of the 
two hyperfine level states, followed by a quantum nonde- 
molition measurement of the population difference of the 
two states, and a sequence specific spin rotations. We 
present the theory of a quantum noise limited interfer- 
ometer, with an atomic cloud placed in one of its arms, as 
a device to be used for the generation of a spin squeezed 
atomic sample. We report the first experiment on non- 
destructive interferometric measurement of the atomic 
population with the sensitivity approaching the quan- 
tum limit. The experimental results are obtained with 
an atomic ensemble in a thermal equilibrium of the two 
hyperfine ground states. Therefore, in the experiment we 
do not measure the projection noise but rather the atomic 
population fluctuations. The goal of the experimental 
part is to show that a nondestructive measurement of 
the atomic population with sensitivity of the order of 
y/Nat, corresponding to the projection noise sensitivity, 
can be achieved with our methods and hence the feasi- 
bility of the method. A nondestructive measurement of 
the atomic level population using phase contrast imaging 
has been reported in 0- 

The paper is organized as follows. In theoretical 
Sec.m we introduce the pseudo-spin in the Bloch sphere 
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picture. We proceed to derive the equations govern- 
ing the interaction of the probe field with the Caesium 
atoms, followed by an analysis of the noise contributions 
to the interferometric signal. The next two theoretical 
sections deal with the effects creating and counteracting 
spin squeezing and we end the theory by illustrating qual- 
itatively how the squeezed pseudo-spin can be incorpo- 
rated into the clock operation with the aim of improving 
its precision. 

In the experimental section we start out with the de- 
scription of our setup and continue to document its op- 
erational properties, with emphasis on the interferome- 
ter noise and methods to supresion unwanted noise. We 
present the results of a measurement of the phase shift 
due to the atoms in Sec. 171 and in Sec. 1111 El we show 
the results of a measurement of the atomic population 
noise of our cold atoms with the sensitivity of \fN at . We 
conclude in Sec. IIVI and present the outlook towards the 
future implementations and improvements of the experi- 
ment. 



to the conservation of the number of atoms, is repre- 
sented as a trajectory of {j) on the surface of a sphere. 
The population difference between the two atomic lev- 
els is then given by the projection of j on the polar axis 
(j z ), whereas the projection onto the sphere's equatorial 
plane gives information on the coherent superposition of 
the two atomic states |3) and |4). 

B. Atomic phase shift 

The initial step is to investigate how an atomic pseudo 
spin can influence the phase of an optical probe field near 
resonance on a transition between the hyperfine ground 
states and an excited state. To this end, we start by writ- 
ing up the complex index of refraction imposed on off- 
resonant light by a sample of cold multilevel atoms. We 
consider the alkali D transition J — > J' between states 
having total electronic angular momenta J and J' . The 
index of refraction is given by |10| 



II. THEORETICAL DESCRIPTION 

A. Pseudo spin in the Bloch sphere representation 

Let us first introduce the two level atom formalism to 
describe the hyperfine ground levels of Caesium. Con- 
sidering a two level atom k with the states |3) and |4) , 
the density matrix elements are — \i){j\k (i,j — 3,4). 
Note, that in atomic clocks the two levels are the m = 
magnetic sublevels but in much of the discussion below 
the magnetic state of the two levels is not important. 
Such a two level atom can be described in terms of a 
spin one half system with the pseudo spin operators de- 
fined by 

Jx = \{P^+Pli) , 

Jy = Y^43"P34) , (1) 
Jz = 2^44 -P3S) , 

where Jx> jy and Jz are the projections of the angular 
momentum operator j k on the x, y and z axes, respec- 
tively. 

For an ensemble of N at atoms, we define the col- 
lective angular momentum operators by j x — j* , 
for the x component, and similarly for the other ones. 
These operators fulfil the angular momentum algebra 
[iijjj] = i £ ijijh where Siji is the Levi-Civita tensor, and 
are useful for illustrating the evolution of the atomic 
quantum state using the Bloch sphere representation. 
This representation is obtained by plotting a vector which 
x, y and z components are given by the mean values of 
j x ,jy and j z , respectively. From this picture, the conser- 
vation of (j) 2 = (jx) 2 + {jy) 2 + {jz) 2 1 which is equivalent 



»A-l=g^(2J+l) (2) 

where I and F are the nuclear and total atomic ground 
state angular momenta respectively, and the primed 
quantum numbers refer to the excited states. We have 
also introduced Np for the atomic density in the level 
with angular momentum F, Appi = wpp - w for the 
detuning of the probe light from the F — > F' transi- 
tion, the atomic linewidth 7 and finally the wavelength 
A, assumed to be common for all transitions making up 
the considered D line. Eq. (J2J is valid for a polarized 
probe interacting with a currently experimentally realis- 
able unpolarized atomic ground state so that the popula- 
tion density in the gound state magnetic sublevel \F,mp) 
is Np^mp = Np/(2F + 1), and we have assumed detun- 
ings small enough to have \App> | <w. For the Caesium 
D2 line (J' = 3/2) of relevance in our experiment, we 
have F = {3,4} and F' = {2,3,4,5}, as shown on Fig. 
|TJi). As we will see, the phase-shift associated with the 
index of refraction @ carries the relevant inforamtion 
about the z-component of the pseudo-spin, and can be 
measured using an interferometer as depicted in Fig. [5] 

Eq. is linked to the pseudo-spin, when we de- 
scribe the population of the two hyperfine ground states 
in terms of the pseudo spin component j z . 

If we consider the situation where both hyperfine 
ground states are close to being equally populated, (N3 = 
N4,), then in the pseudo-spin language we will have 
{j z ) = 0, and let us say only {j x ) = j with a nonzero 
mean value. In the Bloch sphere representation, this sit- 
uation corresponds to a vector in the equatorial plane 
as in Fig. |3Jd). 111 this situation, the atomic variance of 
j z is the same as that of j y and equal to N at /4 [Til . 



3 



a) 



251 MHz 



6 2 P, 



201 
151 " 



'33CS 



.F'=5 

4 

_ 3 
■ 2 



852 nm 



6 2 S« 



9.192 GHz 



. F=4 



b) 



0.15 



0.1 



0.( 



0.C5 



0.1 



0.15 



4260 4280 43JXT 4320 4340 4360 



2000 



4000 



6000 



8000 



AjMHz] 



Figure 1: (Color online) a) Diagram of the Cs hyperfine levels included in the D2 line, b) Theoretically evaluated phase shift 
of the probe as a function of the detuning A45 from the 6S 1 /2{F= 4) — > 6P3/2(F" ' = 5) transition. 



This can be depicted on the Bloch sphere [Fig.[3))] by an 
uncertainty disk at the tip of, and perpendicular to the 
mean value of the angular momentum vector. It is well 
known 0, Q that by performing a QND measurement of 
j z this quantity can acquire a value more well defined 
than that corresponding to an ensemble of independent 
atoms and thus spin squeezing of the pseudo spin vector 
can be achieved. 

In our case, the QND measurement will be performed 
by monitoring the optical phase shift of the off-resonant 
probe interacting with our Cs atoms on the D2 line. This 
phase shift is given by 4>a = fcoffiejnA — 1}, where / is 
the physical length of our atomic sample and ko is the 
optical wavenumber. Using the Eq. we find 
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where 0o = 2 ^ and we have introduced the parame- 
trization N 3 =*jV(l - (3)/2 and iV 4 = N(l + 0)/2, N 
being the overall atomic density and (3 — (N4 — N 3 )/N — 

(jz)/j- 

Using the hyperfine splittings listed in Fig. ^1) and 
inserting the relevant values for the 6J symbols, we 
find by solving J2| a zero phase shift of the probe at 
A /2tt = 4312 MHz relative to the F= 4 -> F'= 5 tran- 
sition. At this detuning the phase shifts from the two 
ground state transitions to the excited state hyperfine 
manifold cancel for equal populations (/3 = 0), as is illus- 
trated in Fig. ^3). Therefore, at Ao any excursions of (3 
will result in an optical phase shift proportional to /3 and 
hence information can be obtained about the collective 
atomic pseudo spin j z and, in particular, the quantum 
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Figure 2: (Color online) The Mach-Zehnder interferometer 
with loss sources and associated input fields indicated. The 
atoms are considered to be in the upper arm. 



fluctuations of this observable 01 . Since the latter are 
the manifestly quantum features of the collective atomic 
pseudo spin observable, a nondestructive measurement 
with the sensitivity at the level of atomic quantum fluc- 
tuations will fix (j z ) to the recorded value at the expense 
of measurement induced back action noise in the orthog- 
onal j y observable. The state determination will, among 
other things, be limited by the accuracy of the measure- 
ment and hence it is important, for the estimation of the 
degree of spin squeezing achievable, to evaluate the lim- 
iting noise sources of our phase measurement. Also rele- 
vant to our study of spin squeezing is the degree to which 
the probe excites transitions in the atomic medium. Ob- 
viously, such excitations will partially cancel the effect 
of the QND measurement and therefore, it may impose 
limitations to the achievable degree of spin squeezing. 
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C. Calculation of the signal-to- noise ratio for the 
quantum noise limited interferometry 

For the monitoring of the atomic phase shift, we con- 
sider the experimental situation illustrated in Fig. [21 
where a Mach-Zehnder interferometer is placed around 
the atomic sample. After the interaction with the atoms 
the transmission of the interferometer arm is 77 and the 
mode overlap at the second beam splitter is y/V. The 
operator corresponding to the input probe field is desig- 
nated by a = al + (x + iy)/2, where a is its real mean 



value and, x and y are the fluctuating quadrature compo- 
nents. The input photon flux is then <£> = a 2 . Similarly, 
we introduce the vacuum fields bk — (xk + iyk)/2 with 
k = 1 to 4, mixing via the loss processes. At the final 
beam splitter we must consider two orthogonal spatial 
modes due to the non perfect mode overlap. Clearly, 
these modes cannot interfere optically but they will how- 
ever add coherently in the detector photocurrents. For 
the photon fluxes impinging on the two detectors we ar- 
rive at 



d\di = ~[1 +77 - 2^/rjV cos(4>)]a 



d\d 2 = - [l+r] + 2\Jr]V cos((^)]a t a 
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+ h.c. , (5) 



r 



where the phase difference (j> — 2Tr( y J re ^n(L)dL — 
J be n(L)dL^ I A is 27r/A times the difference of the in- 
terferometer arms' optical path length, i.e., the integral 
of the index of refraction over the respective arm. From 
this expression of <f> it is clear that the phase can shift 
because of a change in either the path lengths or the in- 
dex of refraction in one of the arms, or because of a shift 
of the wavelength of the probe. Moreover, in the probe 
arm the atoms can change the index of refraction and 
induce a phaseshift a , so that we write = + a for 
the total phase shift. The last contribution provides the 
monitoring of the atoms, while the first three add noise 
to the measurement. 

Now, we find the mean photocurrent difference to be 



-) = {d\di)- 



d\d 2 



■qV cos{(t>) , (6) 



in units of elementary charge. The visibility of our inter- 
ference fringe is found from the single detector photocur- 
rent and is given by 
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In the symmetric case, where 77 = 1, this reduces to V = 
VV as expected. 

The fluctuating part of i_ is now calculated by lin- 
earizing around the mean value, Si- = i- — (?_), and 
remembering that only a) a has nonzero mean. As a re- 
sult we find 



r)V[cos(4>)x + sin(0)xi 



V(l - 77) 



cos(0)i 2 + sin(0)y 2 ] 



-V3 



1 - V 



cos(0):E4 + sin(0)j/4] 



(8) 



All the field operators in (JSJ are uncorrelated and con- tribute with 2B, where B is the bandwidth of our mea- 
sequently for a coherent state input all operators con- 
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surement ^ij. From this, we find (Si-) 2 oh = Ba 2 (l + T]), 
which is just 2B times the total photon flux, $, trans- 
mitted through the interferometer. 

Let's assume that there are no atoms in the probe arm 
so that 4>a = 0. To be sensitive to small phase shifts, we 
use a second laser far away from the atomic resonance 
to lock the interferometer at the side of the interference 
fringe. With this procedure applied to the system, we 
set the residual phase <f> equal to 7r(l/2 + to) with to = 
0,±1, which has the following consequences: Even if 
our input state is not coherent, or in other words, we 
are probing the atoms using a noisy laser, we will find 
the amplitude noise of the probe laser [the first term in 
©] to be considerably suppressed due to the balanced 
detection. However, the laser phase noise will remain 
important. We model this noise as an excess noise of 
the vacuum inputs x\ and §2, interfering with the probe 
and contributing with a variance (1 + TV) relative to the 
vacuum state noise while 2/3, £4 and ^4 still are at the 
vacuum noise level. 

Considering now the presence of atoms, their contri- 
bution to the phase noise is denoted (50a) 2 , which like 
the laser phase noise is normalized to the probe vacuum 
noise level. Incorporating these values into JSJ) and tak- 
ing into account the transmission of the interferometer, 
we arrive at 



(Si. 



2B + V& tV+( — (5<f> A ) 2 cos 2 (0 A ) , 



1 



n 



again in units of quantum noise of the transmitted probe 
and for the case 4> = 7r (l/2 + to) + ^a- 

The excess noise TV can be suppressed by operating the 
interferometer in the white light position since 



TV = 



Suj 2 



(Slo 2 



-(fc Ai) 2 



(10) 



where AL is the optical path difference between the two 
arms of the interferometer, Suj 2 is the laser frequency 
noise and (Sui 2 ) q — uJq/o 2 is the quantum level of fre- 
quency noise. 

We can also straightforward generalize Eq. @ to de- 
tectors with less than unity quantum efficiency, e. In this 
case we get 



(Si-)" 



= 2B + eV<& 



n 



l + r) 



cos 2 (</>a) 



(11) 

Having addressed the optical noise contributions we will 
now turn to consider the atomic imprint on the probe 
phase noise. 



this is clear because a higher signal to noise yields more 
knowledge of the atomic spin observable, hence it be- 
comes better defined and higher degree of squeezing of 
that observable is achieved. 

Assuming that we have a white light interferometer, 
the noise contributions of relevance here are the quan- 
tum fluctuations (Si) 2 of the phase of the probe pulse 
and the electronic noise (Si) 2 of our photodetectors. The 
probe pulse is characterized by the duration r and the 
photon number $r. Since a highly coherent laser is be- 
ing used to generate the pulse, we assume it to be Fourier 
limited [lj|, 2ttBt = 1, and then 



(12) 



where we have used the photon flux detected from the 
probe arm r]<&/2 as reference for the shot noise level. 

The electronic noise can be described by the noise 
equivalent power (NEP), P e , so that 



(Si) 2 e = (P e \/hc) 2 /2irT 



(13) 



h being Planck's constant. 

If we ignore the electronic noise [(Si) 2 , <C (Si) 2 ,] that 
is only relevant when feedback schemes are involved, we 
find that the signal to noise ratio of the measurement is 
given by 



eV 



2r, 



(l + r?) : 



-7t$t(50a) 2 cos 2 (0a) 



(14) 



which is related to the degree of spin squeezing ^5| , as 
described below. 

The atomic contribution to the phase noise is com- 
puted from Eq. ©. We assume that the atoms ini- 
tially are prepared in a coherent spin state for which 
( S P)loh = ( S jl)coh/j 2 = N^ 1 , where N at is the num- 
ber of atoms within the probe volume. Hence, we find 
that 



> 2 P(A) V 
AttA ) 



N at 



(15) 



where A is the probe beam cross sectional area found 
from the probe beam waist wq as ttwq/2 and we have 
defined the detuning function 




A 2 4- -v 2 



(16) 



D. Spin squeezing 



Finally, we find the ratio of the phase noise from the 
atoms to the quantum phase noise to be 



In the context of spin squeezing a high ratio of atomic 
spin noise to optical quantum noise is desired. Intuitively 
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For our pulsed measurement we integrate i_ over the 
pulse duration and analyze the statistics of collections 
of pulses. This sets an upper limit to the frequency of 
the fluctuations, that can be observed, at approximately 
t _1 . The lower limit is simply set by the time over which 
we collect the integrated pulses. Since the atomic noise 
spectrum is not white, we stress that (|17H is only valid 
in as much as we match our pulse spectrum to cover the 
atomic noise spectrum. 

With respect to the spin squeezing, we will be using 
the following definition taken from 



t = (Sp) 2 N at = (Sf z )N at /f 



(18) 



measurement of the first pulse 01 • If we have created a 
spin squeezed ensemble it will reveal itself via the reduced 



variance 



S(a\ — CL2) 2 < Sa 2 + 8a\ — 28a 2 = 25a 2 , 



(21) 



The quantum nature of the atomic fluctuations can be 
verified by showing that 8a\ and Sa 2 are equal and grow 
linearly with N at . From this experimental considaration 
it is already clear that the sensitivity of the detection 
apparatus must be large enough for the N at noise to be 
detected. 



where £ = 1 for a coherent state, £ < 1 for a squeezed 
state and £ = 00 for a thermal state. This is not the 
only way to define the spin squeezing parameter, but 
this definition characterizes the quality of the state in 
a spectroscopic measurement, i.e., it is a measure for the 
increased sensitivity to rotation in the squeezed direction 
on the Bloch sphere |l6|. 

The degree of spin squeezing can be shown 

dEl to 

be related to n through 



1 



1 



(19) 



and thus, the squeezing imprinted by the measurement 
onto the z-component of the spin is 



E. The relation between the flgure-of-merit and 
atomic decoherence 



The effect counteracting the spin squeezing is the inco- 
herent transfer of atoms from the spin squeezed state via 
optical excitation and spontaneous emission to a mixed 
ground state. This excitation happens with a pulse inte- 
grated rate p e 



Pe 



er(A)$7 



(22) 



where the absorption cross section for the probe is 
er(A) = (A 2 /37r)£(A) with the linewidth function given 
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) sq 



1 N n 



z)coh 



41 + K 2 



(20) 



so that, in order to perform a good QND measurement 
and hence, to achieve a high degree of spin squeezing, 
we must have k large compared to unity. It is natural 
therefore to call k the figure of merit of the QND inter- 
action. Below we show how k can be expressed via easily 
accessible experimental parameters. It is important to 
note that as long as the spontaneous emission rate over a 
pulse is negligible the minimal uncertaintystate will be 
preserved by the phase shift measurement pi EH 13' and 
in this case the variance of the conjugate spin component 
will become (Sj y ) 2 q = jN at (l + k 2 ). 

So far, we have ignored the electronic noise in the above 
calculations. It is however important if feedback schemes 
should be app lied in order to either enhance the spin 
squeezing 20], or if we wish, to rotate the mean spin 
direction according to our measurement with the goal 
of obtaining a specific spin squeezed state. The latter 
is relevant if the spin state should be employed in, e.g., 
atomic clocks, where the a Bloch vector in the equatorial 
plane is desired plj . 

The degree of spin squeezing can be measured exper- 
imentally by sending pairs of probe pulses through the 
atomic sample, integrating these pulses and storing the 
resulting areas, a± and ai- The variances Sa 2 = 8a 2 , will 
set the level of the atomic quantum noise, while the vari- 
ance of the pulse difference 8[a\ — Q.2) 2 will yield informa- 
tion of interatomic correlations created by the quantum 
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(23) 



The above cross section assumes both ground states 
having equal populations and we find that the pulse in- 
tegrated excitation rate is related to K through 



= eV 



V 



X 2 V(A) 2 



(1 + V ) 2 4A £(A) 



COS 2 (<j) a) NatPe 



(24) 



To highlight the relevant physical parameters this 
equation can be cast into a compact form. If we assume 
the visibility V, quantum efficiency e as well as the trans- 
mission of the interferometer 77 all equal to 1, then in the 
limit of large detunings A 3> 7 we have 2?(A) 2 /£(A) w 1 
and the DC-phase shift becomes negligible so that also 
cos(</>a) ~ 1, and Eq. simplifies to 



(25) 



where we have introduced ag, the atomic optical density 
on resonance. It is clear that, in order to achieve strong 
spin squeezing we need a large ao since we wish to keep 
p e small to maintain the nondemolishing character of the 
measurement. 
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Following the excitation the atoms decay sponta- 
neously to a spherical spin state, i.e., a state character- 
ized by zero expectation value of all spin components 
(ji) = (i = x,y,z). This dissipative evolution of the 
atomic state leads to the additional noise into the final 
state, resulting in less efficient spin squeezing. For small 
p e the degradation is of the order of p e . A detailed analy- 
sis of the spin state resulting from the QND measurement 
and a certain amount of excitation is presented in |22j | . 



F. Spin squeezing in clock operation 

In the preceeding sections we have shown how a non- 
destructive optical phase shift measurement reduces the 
noise of the atomic pseudo-spin z-component. We will 
now put this QND measurement into the context of the 
clock operation, which is discussed qualitatively in the 
Bloch sphere representation. More detailed accounts of 
the atomic clock operation in this picture can be found 
in 0. 

The full protocol including the standard clock sequence 
can be viewed as follows: Initially we use optical pump- 
ing to prepare the atoms in a coherent spin state with 
(j z ) = -N a t/2 and (j x ) = {%) = 0, as shown in Fig.^Jri. 
This is the situation where the atomic variances of j x 
and j y are both equal to N a t /4 as depicted by the uncer- 
tainty disc on the Bloch sphere [Fig.^)]. Next step is to 
apply a classical ir/2 pulse using an RF-magnetic field, 
which corresponds to the first it/2 pulse in the Ramsey 
spectroscopy sequence. This pulse brings the angular mo- 
mentum vector to the equatorial plane as illustrated on 
Fig. |3Jd). For a standard atomic clock the spin would 
be allowed to precess in the equatorial plane of the Bloch 
sphere [Fig.Efc)] until the second ir/2 pulse in the Ramsey 
sequence is applied [Fig. |3Ji)] and the atomic population 
difference is detected. Instead of this, we proceed from 
the state depicted on Fig. ISJd) to perform a QND mea- 
surement of the population difference j z using an optical 
field. As argued above, this measurement reduces the 
uncertainty of the operator j z at the expense of an in- 
creased uncertainty of j y while preserving the minimal 
uncertainty state and so, the atomic sample is prepared 
in a spin squeezed state. This is the state illustrated in 

Fig. Eft). 

The degradation of the spin squeezing that can be 
caused by spontaneous emission will cause an increase 
in the size of the uncertainty disc, as well as reduce the 
length of the Bloch vector leading to the loss of contrast 
in the clock signal. However, both effects are of the or- 
der of the pulse integrated rate of spontaneous emission 
which can be kept small, simply by letting the probe be 
far detuned. 

Additionally, we must update the atomic state based 
on the result of the QND measurement j^jj. Depending 
on the outcome of the measurement the mean value of 
the pseudo spin vector will be shifted away from (j z ) = 0. 
This deviation will be corrected for by application of a 
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Figure 3: (Color online) The Bloch-Sphere scheme for the 
state preparation (a,b) QND measurement (e,f,g) and the 
Ramsey spectroscopy for a coherent state (c,d) and a spin 
squeezed state (h,i), with the corresponding simulations of 
the projection noise on the Ramsey fringe, (j) and (k) respec- 
tively. 



short RF-pulse that will shift the vector back into the 
equatorial plane of the Bloch sphere [Fig. Of)] ■ 

Since for the Caesium clock it is important to have 
reduced noise during the precession of the phase com- 
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ponent j y in the equatorial plane, we rotate the pseudo 
spin vector around the x axis with a ir/2 pulse [Fig. [3k)], 
which effectively interchanges the j y and j z components. 
For the 7r/2 rotation to revolve around the x-axis, the 
RF-magnetic field, applied in this step, must be phase- 
shifted by 7r/2 w.r.t. that the of the first Ramsey 
pulse [Fig. Ob)]- The following steps, represented by 
Figs. and 0), correspond to the standard method 
for Ramsey spectroscopy in the clock operation, but now 
using a spin squeezed state. 

If we compare Figs. [31) and[3Jt) we clearly see that we 
would gain in signal to noise ratio of the spectroscopy 
signal from the projection noise limited measurement of 
the population difference (j z ) in the clock transition [l6| . 
performed at (lo — ll>q)T — ±7r/2. 



III. ATOMIC NOISE MEASUREMENTS 



Probe Laser 



BS 

50/50 FPD 



SA 



1MHz 



Reference Laser DBM 

Feedback to Amp 
Probe Laser 
Controller *~ 



Figure 4: (Color online) Experimental setup employed to lock 
the probe beam. The elements included in the sketch are: 
BS - 50/50 beam splitter; FPD - fast photodetector; SA - 
spectrum analyser; DBM - double balanced mixer; LPF - low 
pass filter; Amp - amplifier. 




In this part we aim at showing that our apparatus has 
the sufficient sensitivity to measure the atomic noise and 
along the way we analyze the various considerations and 
precautions necessary for reaching the goal. 

A. Magneto-optical trap 

The atomic sample is prepared in a standard six beam 
Cs magneto-optical trap (MOT). We are able to trap 
around 3 x 10 s atoms, when loading the MOT from back- 
ground Cs vapor with sufficiently high partial pressure 
(around 10~ 7 mbar). The red detuning of the trapping 
laser is set to 15 MHz. The cloud volume is approxi- 
mately 6 x 10~ 3 cm 3 , which at best can yield a resonant 
optical density of 13. Due to high background pressure 
the lifetime of the trap is only around 20 ms. This short 
lifetime is convenient for acquiring statistical data be- 
cause it allows for quick refreshing of the atomic sample. 

B. Frequency locking of the probe laser 

The probe laser is locked and blue detuned from the 
atomic transition &Si/ 2 {F = 4) — > 6P 3 / 2 (F' — 5) by a 
specific detuning A variable from a few MHz to a few 
GHz. The experimental setup used to lock the probe 
laser in this way is shown in Fig. 

Two lasers with a specific relative frequency separa- 
tion can be locked in a number of ways [23L |24| . We 
use a reference laser locked to the atomic transition 
6Sy 2 (F = 4) -> 6P 3/2 (F' = 5) by FM saturation spec- 
troscopy. The reference and a fraction of the probe laser 
beams are mixed at the beam splitter BS and their beat 
note is measured using a Newport fast photodetector 
(model 1480) with 15 GHz bandwidth. The produced 
RF signal has a frequency component corresponding to 
the probe detuning from the above specified atomic tran- 
sition. This signal is monitored with the spectrum an- 



alyzer SA, as shown in the figure. From the RF signal 
obtained in such a way, we can generate an error signal 
to lock the probe laser with the desired detuning from 
the atomic transition. To accomplish that, the current 
of the probe laser is FM modulated with a modulation 
depth of 1% using a 1 MHz sinusoidal waveform that is 
also utilized as a local oscillator for the double balanced 
mixer DBM in Fig. 0] The output from the mixer is a 
DC signal that is low pass filtered and amplified before 
feeding it back to the probe laser controller. 



C. White light fiber interferometer 

The interferometer as shown on Fig. [SI is a Mach- 
Zehnder type made of single mode optical fibers. The 
motivation for using fibers instead of free space prop- 
agating beams is the enhanced mechanical stability as 
well as excellent mode overlap of the interfering beams 
in single mode fibers. The field in the input fiber enters 
a 50/50 coupler CI and is split into a reference arm and 
a probe arm surrounding the atoms to be probed. The 
field in the probe arm exits the fiber and with the lens 
LI it is focused at the center of the MOT with a beam 
waist of 20 /mi. After passing through the atoms, it again 
enters the fiber and is combined with the field from the 
reference arm at the second 50/50 coupler C2. 

Since we use a non polarization maintaining fiber, the 
field polarization can evolve differently in the two arms. 
Thus, in order to achieve maximal interference visibility, 
we include polarization controllers PC to match the field 
polarization from the two arms at the second coupler C2. 
We note that the coupling efficiency r\ through the air- 
gap containing the MOT is 30 %. With the fiber couplers 
providing a nearly perfect mode overlap of the probe and 
the reference field, i.e., y/V ~ 1, the visibility becomes 
V = 85%. 

The pulsed probe signal is detected with a balanced 
detection scheme [2fJ using the low noise photodiodes Dl 
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Figure 5: (Color online) Sketch of the setup of the interferometer with following elements: BS - 50/50 beam-splitter; CI & C2 
- 50/50 fiber couplers; PCI, PC2 & PC3 - fiber polarization controllers; LI & L2 - achromatic lenses; Fl & F2 - interference 
filters transmitting @ 852 nm; Dl & D2 - Hamamatsu low noise, high gain photodiodes; D3 & D4 - photodetectors; and several 
half wave plates A/2 and collimating lenses for fiber coupling, ih is the locking signal, whereas i_ = i\ — i-z is the probe signal. 



and D2. From the integral of the photocurrent z_ over 
the pulse duration, we extract the area which corresponds 
to the difference of the signals from the two arms. The 
mean value of the difference gives the DC-phase shift 
4>, and the variance gives information about the phase 
fluctuations. 



1. Locking the Interferometer 



2. Shot noise limited interferometer 

With all the measures for removing undesirable noise 
contributions implemented, we now try to gauge their 
effectivenes. To get a measure of the instability or noise 
of the interferometer signal we measure the variation of 
the area from one pulse (a^) to another (aj+i). This is 
done by determining the two point variance |27| cr 2 (ro) 
which can be defined to be 



To reduce thermal and acoustic drifts of the interfer- 
ometer, we lock it by means of an off-resonant CW laser 
that propagates through the interferometer simultane- 
ously with, and in the same direction as the probe beam. 
The locking beam is several nm away from the atomic res- 
onance and therefore it is not affected by the cold atoms. 
At the output the locking beam and the probe beam are 
separated by the interference filters Fl and F2. The bal- 
anced detectors D3 and DA provide an error signal which 
controls the piezo adjusting the length of the probe arm. 

In order to cancel the amplitude noise term in the 
Eq. ©, the interferometer needs to be locked so that 
when the cold atoms are absent (0a = 0) the interference 
signal for the probe is at half fringe, i.e., cf) = 7t(1/2 + to) 
with m = 0, ±1, .... However in this position the inter- 
ference signal is also most sensitive to the phase noise 
of lasers. We use semiconductor lasers that are charac- 
terized by strong phase fluctuations with a wide band of 
frequencies for both probing and locking 26]. We have 
measured their linewidth to be approximately 500kHz. 
As discussed in Sec. Ill CI in order to suppress the ef- 
fect of the phase noise we lock the interferometer at the 
white light position, corresponding to a nearly zero path 
length difference |2g. We use a regular broadband LED 
as a white light source to determine roughly the white 
light position. 



^ M-l 

a2(To) = 2 (Af-l) ^ ~ ^ ' (26) 

where tq is the temporal pulse separation and M is 
the number of pulses in our measurement. The pulse 
sequences used to compute c 2 (to) are composed of sev- 
eral thousand pulses each of 2 /is duration. From these 
sequences we can extract the two point variance on 
timescales comparable with the pulse duration and up 
to two orders of magnitude in r. The results correspond- 
ing to this measurement are shown in Fig. [5J 

If the interferometer noise would be purely white, the 
two point variance would stay constant on all timescales. 
Naturally, temperature drifts would cause the variance 
to rise on larger timescales than those we measure. We 
observe that the two point variance fluctuates with a pe- 
riod of 200 ms corresponding to 50 Hz line noise. On 
the scale of Fig. this is seen as a slow rise of the 
curves. However, on the /is timescale this line noise is 
of no importance. One interesting feature is seen from 
the upper trace of Fig. [fj] By modulating the laser diode 
current and thereby, the frequency of the probe laser, 
we see that cr 2 (r ) oscillates with a period corresponding 
to the modulation frequency. This in turn means that 
the interferometer is sensitive to frequency changes and 
therefore, isn't exactly in the white light position. By 
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Figure 6: (Color online) The two point variance extracted 
from measurements done with pulse separation to = 20 [is. 
Lower trace: Interferometer in white light position with probe 
laser locked and the variance increases on larger timescales 
due to the 50 Hz line noise. Upper trace: Interferometer out 
of white light position with probe laser frequency modulated 
at 5 kHz, which is directly reflected in the variance by a 200 /is 
oscillation period. Additional phase noise from the laser raises 

the level of the minima ( ) with respect to the lower trace, 

because the laser is not phase locked and the interferometer 
is not in the white light position. 




100 1000 

Probe power [nW] 



Figure 7: (Color online) Noise in the amplitude x (o) and 
phase y (•) quadratures of the probe light. Fits to the data 

in log-log scale give slopes of 1.2 ±0.2 ( ) for the amplitude 

and 1.2 ± 0.4 ( ) for the phase quadrature. 



adjusting the path length difference of the interferometer 
arms it is possible to come in to the white light position 
and consequently, the oscillation disappears from the two 
point variance as seen on the lower trace of Fig. This 
method ensures a fine white light alignment to within 10 
jUm, so that we can have m < 2 and fulfill the condition 
<j> < 5tt/2. Since it is of less concern what is the exact 
value of m, for the sake of simplicity we assume it to be 
zero. 

It is crucial to find how the interferometer noise, as ex- 
pressed by the variance of the probe pulses, depends on 



the power of the probe laser. If the interferometer is shot 
noise limited, then according to Eq. (|llf) we would expect 
the noise to scale linearly with the probe power. Apply- 
ing the same pulse sequence as described above, we first 
try to determine the noise in the amplitude quadrature. 
For that purpose, we block the probe arm of the inter- 
ferometer and just observe the noise of the transmitted 
beam. In Fig. a fit to the data shows that the am- 
plitude noise depends linearly on the probe power. This 
shows that our detection is shot noise limited and in ad- 
dition to that, we have found the shot noise level for later 
reference. 

Next, we unblock the probe arm and look at the phase 
noise of the interferometer, which in the white light po- 
sition should be dominated by shot noise. This is also 
shown in Fig. [7] and one sees that the dependence on the 
probe power is linear to within the uncertainty of the 
fit. On this, we conclude that the interferometer is shot 
noise limited in the probe power range 0.05 /xW - 2 /iW. 
Moreover, Fig. [7| shows that the quantum noise is simi- 
lar in the two quadratures, as one would expect from a 
coherent probe field. 

We are now in a position to apply the interferometric 
detection to a sample of cold atoms. As a start we show 
how the interferometer can be used to nondestructively 
measure the number density of the trapped atomic cloud. 



D. Interferometry with cold atoms 

We begin with balancing the interferometer at half 
fringe in the absence of trapped atoms. Trapped atoms 
cause a DC-phase shift of the probe pulse. The mea- 
sured signal must then be corrected for the absorption of 
the probe. From the corrected signal inc, the DC-phase 
shift </>a is deduced as 4>a = arccos (idc)- Measuring </>a 
as a function of the detuning, in the vicinity of F = 4 
level, we observe the dispersive behavior shown in Fig. [S] 
Fitting the experimental data points using the Eq. 
we can determine the density of atoms 
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4.3 10 9 cm~ 



(27) 



where C is the parameter given by the fit and is propor- 
tional to 4>q. The estimated number of probed atoms is 
then 5500. 

Notice that the theoretical fit on Fig. [H] overestimates 
the amplitude of the 6S 1/2 (F = 4) 6P 3/2 (F' = 3,4) 
transitions. This can be explained in the following way. 
We have accounted for the reduction in photon flux $ 
and thereby the amplitude of the interference fringes 
as a result of absorption. However, we did not con- 
sider the optical pumping of the atoms due to the ab- 
sorbed light, and the consequent reduction of the respec- 
tive atomic density. The reduction of the atomic den- 
sity at each resonance will depend mainly on the res- 
onant term of the lineshape function T>(A), where the 
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Figure 8: (Color online) DC-phase shift due to dispersion by 

the atomic cloud with fit ( ) to Eq. A closer fit ( ) 

is obtained by letting the line strengths be independent, and 
thus not given only by the value of the 6J-symbols. 



weight factor of the resonant terms should give the re- 
lationship between the absorption probabilities at the 
three resonances. However, for the QSi/ 2 {F = 4) — > 
6Ps/2(F' — 3,4) transitions we are pumping atoms over 
into the F = 3 ground state, making them insensitive to 
the probe at the current detuning. On the other hand, 
the 6Si/ 2 (F= 4) — * 6P 3 / 2 (F' = 5) transition is a cycling 
transition since the decay to the F — 3 ground state is 
not allowed. Therefore, we do not depump the F = 4 
ground state at this resonance. Consequently, the fit 
to Eq. (|3J) should avoid points close to the depumped 
QS 1/2 (F = 4) -> 6P 3/2 (F' = 3,4) resonances, and the 
number C used in Eq. (|27[) is indeed obtained from such 
a fit. 

Closer theoretical calculations that also take stimu- 
lated emission and the branching ratios of the excited 
level decay into account, allow us to determine time evo- 
lution of the level populations for a given detuning and 
power of the probe. From this the average number of 
atoms that during the pulse have been pumped over into 
the F = 3 ground state can be calculated, and used 
to estimate the reduction of the DC phase shift due to 
depumping. With this compensation applied to the data, 
the ratios between the amplitudes of the DC phase shift 
at the three resonances, come very close to values pre- 
dicted by the transition strengths. 

We conclude that in the vicinity of the F — 4 — > F' = 5 
transition the probe does not perturb the atomic popu- 
lation with our measurement conditions. 



E. Poisson noise of atomic population 

As we can see from Eq. (|2Jl, at constant detuning the 
phase shift is proportional to the density of atoms N. 
Therefore, given a fixed detuning we can use the DC- 
phase shift as a measure for the number of atoms probed 
in the MOT [H |29|. We vary the number of atoms 



trapped by varying the background Caesium pressure in 
the chamber. 

The precision of our setup that we have so far reached 
does not enable us to measure at the detuning of 
A45/27T = 4312MHz where the atomic phaseshift can- 
cels (cj)A = 0) for equal populations in the two hyperfine 
ground states. With this in mind, the optimal choice of 
the probe detuning is dictated by the balance between 
the strength of the QND interaction and the strength of 
decoherence processes. As in Sec. Ill Dl the signal refers 
to the photocurrent variance due to atoms and the noise 
is the photocurrent variance due to other noise sources. 
Eq. (fT7|) tells us that the signal to noise ratio k 2 goes as 
the square of the detuning function T>(A) and linearly 
with the probe flux From this relationship, we in- 
fer that a better signal is obtained at small detunings 
from the transitions where 2?(A) is large. On the other 
hand, at small detunings the photon flux would decrease 
strongly due to absorption by the atoms. Moreover, the 
excitations can destroy the atomic coherences. 

For the data presented here, the atomic population 
noise measurements were performed at the blue detun- 
ing of 15 MHz from the F = 4 — > F' = 5 transition. As 
a consequence we measure the noise of the atomic pop- 
ulation in the F = 4 ground state rather than the noise 
of the pseudospin component j z . Using 2 fis long pulses 
of 0.6 /zW power, this yields a pulse integrated rate of 
atomic transitions of p e = 15. Obviously, this is far from 
being a QND measurement and the atomic coherence 
would have been completely lost under such measure- 
ment conditions. Normally, this would also mean that 
many atoms are transfered to the F = 3 ground state, 
but since we are close to a cycling transition, virtually 
all of the atoms remain in the F = 4 ground state after 
the measurement. Therefore, we have effectively relaxed 
the demand for the measurement to conserve the popula- 
tions only, but not the coherences. However, the nature 
of the noise is unaffected by the measurement being or 
not being QND, and in that light the ability to measure 
the noise is a relevant indicator for the feasability of this 
procedure, in spite of the large amount of real transitions. 

Note that the photon flux is limited only by the detec- 
tion, which is shot noise limited up to around 1 fiW. 
With this photon flux of 5 x 10 6 , every atom per- 
forms about 50 absorption cycles with the probe on res- 
onance, whereby their recoil velocity reaches approxi- 
mately 1.5 cm/s, which doesn't pose any problems on 
the timescales of our experiment. 

For this measurement we employ a slightly more com- 
plex scheme of pulses. Each sequence consists of 3 pulses 
of a duration of 2 /is and equally separated by 10 ms. All 
three light pulses have the same 0.6 /iW incident power 
that lies in the range where our detection is shot noise 
limited. Before each sequence we load the MOT and 
turn off the magnetic field just before the first of the 
three pulses. This means that the first pulse 1 probes 
the atomic cloud. With no trapping force to contain it, 
the MOT has long decayed at the arrival of the second 
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Figure 9: (Color online) Phase noise induced in the probe 
light from the interaction with cold atoms. The DC-phase 
shift is used as a measure for the number of atoms. 

pulse 2. In this way, the pulses 2 and 3 probe the chamber 
without any atomic cloud and are used as a reference. We 
repeat this sequence several thousand times to obtain the 
statistical information on the pulse areas a±, a 2 and 0,3. 
We subtract the pulse areas pairwise g?i_2 = cti — &2 and 
c?2-3 = &2 — 0-3, and calculate the mean values di-2, ^2-3 
and variances (<5e?i_2) 2 , (^2-3) 2 , over all the recorded 
sequences. 

Let us address the significance of these variances. For 
the first pulse the atoms are trapped in the MOT and 
thus, the pulse area variance (i5ai) 2 is given by Eq. i|ll|) 
where, as mentioned earlier, we may neglect the classical 
phase noise AT due to the white light alignment of the 
interferometer 

(Sa,) 2 = (5i) 2 p + e 2 Vt> 2 (t^) Wa) 2 cos 2 (0 A ) • 

For atoms in a coherent superposition of two ground hy- 
pcrfine state probed by light tuned in between the two 
states the atomic contribution to the phase noise (5</>a) 2 
is induced by the quantum projection noise as in Eq. (|15|l . 
For atoms in a thermal state probed by light tuned close 
to one of the hyperfine transitions, as in our present 
setup, the atomic contribution to the phase noise is dom- 
inated by the population noise. The population noise 
arises from the statistical nature of the trapping process 
and from the motion of atoms in and out the probe re- 
gion. As we demonstrate experimentally, it is character- 
ized by Poission statistics with the width y/N a t, same as 
for the projection noise. 

For the last two pulses, recorded without any cold 
atoms (0a = 0), there is no noise term from the atoms 
and consequently, the variance of the two areas may be 
written as 

(5a,) 2 = {5i)l , j — 2, 3. (29) 
We see that one ought to be able to extract the shot 



noise and atomic noise from the appropriate pulse ar- 
eas. What hinders it are slow thermal drifts of the inter- 
ferometer that would dominate the variance if not cor- 
rected for. This is why we subtract the pulses pairwise 
to cancel the effect of the thermal drifts. Since the shot 
noise and atomic noise are uncorrelated it is easy to ver- 
ify that for the subtracted areas di_ 2 and d 2 s we have 
(8d^ 2 ) 2 = (5 ai ) 2 + (Sa 2 ) 2 and (<M 2 _ 3 ) 2 = (8a 2 ) 2 + (8a 3 ) 2 
so that we can determine the shot noise and atomic noise 
contributions as 




respectively. 

In Fig. El we show the atomic contribution to the phase 
noise (£0a) 2 as a function of the DC-phase shift <^a, 
which, in turn, is proportional to the number of atoms. 
The linear fit, within the uncertainty, shows that the vari- 
ance of the atomic fluctuations scale as N at . Thus we 
conclude that the white light interferometry is capable of 
achieving the sensitivity at the level of projection noise 
fluctuations. 

IV. SUMMARY AND OUTLOOK 

In this paper we propose a sequence of QND measure- 
ments and spin rotations which allows to circumvent the 
projection noise limit of accuracy of the Cs atom clock. 
The realization of this protocol requires interferometric 
QND measurement of the atomic population with the 
sensitivity at the projection noise level. Towards this 
goal we demonstrate experimentally that the shot noise 
limited fiber optical interferometer at a white light set- 
ting can reach the sensitivity sufficient to detect the pro- 
jection noise under conditions not far from those of the 
QND measurement. 

The main parameter that can significantly improve the 
QND figure of merit, k 2 , is the resonant optical density 
= , which is directly proportional to the number 
of atoms in the probing region. Increasing the optical 
density should allow to reduce the photon number and 
increase the detuning of the probe. Both those measures 
will improve the spin squeezed state preparation. The 
atomic density can be greatly increased compared to our 
present level by applying a far off resonant trap (FORT) , 
where the dipole force from a focused laser beam traps 
the cold atoms |30j . This naturally suggests that the 
QND scheme could be used in an optical lattice clock 
where the number of atoms is preserved and it is possible 
to have high densities. Given the density limitation on 
the clock, an optimal density and atomic number should 
be sought for, in order to improve the clock performance. 

Another improvement lies in the optimization of the 
interferometer, with the aim to reduce losses. This may 
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mean abandoning the fiber interferometer thereby prac- 
tically eliminating coupling losses (77 = 1). There is an- 
other reason why the increased flexibility of a free space 
interferometer can be an advantage. If the atoms are con- 
fined in a dipole trap, the sample will be needle shaped, 
which will have a lensing effect on the tightly focused 
probe beam. In this case it is desirable to be able to in- 
clude compensating lenses into the probe arm. The effect 
of diffraction by atomic samples with different densities 
and geometries on the spin squeezing is analyzed in detail 
in I31I. 
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